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Abstract 

A solution of the nonlinear Klein-Gordon equation perturbed by a 
parametric driver is studied. The frequency of the parametric pertur- 
bation varies slowly and passes through a resonant value. It yields a 
change in a solution. We obtain a connection formula for the asymp- 
totic solution before and after the resonance. 



1 Introduction 

This work is devoted to the problem on a control of a nearly monochromatic 
weakly nonlinear dispersive wave with small amplitude in a strong nonlin- 
ear media. It is well-known that packets of nearly monochromatic waves 
propagate without changing of their shape when the envelope function of the 
packet is a soliton of the Nonlinear Schrodinger equation(NLSE). The soli- 
tary packets of waves would be more suitable for communication in optical 
fibers on a large distance if one can control the parameters of the envelope 
function for such packets. The wave packets with a soliton-like shape form 
for sufficiently large range of initial data. The control of the shape is possible 
by a parametric perturbation of the system. 

Here we propose a new approach for controlling of parameters of the 
solitary packets. In our approach the wave packets are controlled due to a 
slowly passage through a local parametric resonance. 
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In this paper we study the parametrically driven nonhnear Klein-Gordon 
equation. The frequency of the driver varies slowly. It is well-known that an 
envelope function of weak nonlinear wave is a solution of the NLSE [I]- [3]. 
In general case this is valid proposition for the parametrically driven Klein- 
Gordon equation. But the small driving force plays a central role in narrow 
layers where the frequency passes through the resonant value. As a result 
fast and slow variables appear which define the behaviors of the solution. In 
these layers the envelope function of a primary parametric resonance equation 
instead of the NLSE. In the slow variables the passage through the resonant 
layer looks as a jump of the envelope function. Our goal is a connection 
formula for this jump. 

The change of the solution in the resonant layers for a local resonance (not 
parametric) was studied in [1]-[H]- The detail description of the phenomenon 
was presented in [9] . In that case the solution in the resonant layer is defined 
by the Fresnel integral. For parametric resonance the solution was studied 
in [To] for the Mathieu equation. They shown that the solution defines by 
parabolic cylinder functions. In our work we found that the parabolic cylinder 
equation also allows one to obtain the connection formulas for the small 
solutions of nonlinear Klein-Gordon equation. 

In this work we use the singular perturbation theory and matching of 
different asymptotic expansions to obtain the connection formula. Our ap- 
proach is based on the matching method of asymptotic expansions [TT] . 

The structure of the paper is as follows. In Section 2 we formulate the 
result. The numerical simulations are inserted in Section 3. Section 4 con- 
tains the formal construction for the asymptotic solution out of the resonant 
layers. In Section 5 the asymptotic solution is constructed in the resonant 
layer. In Section 6 we match the main terms of the asymptotic expansions. 

2 Statement of the problem and main result 

We study the Klein-Gordon equation with a cubic nonlinearity 

d^U-dlU+(l+efcos(^^^^^^^^y + ^U' = 0, 0<£«1. (1) 

Here 7 and / are constants. The phase function S{y, z) and all derivatives 
of S{y, z) are bounded. 

Our goal is to obtain an asymptotic solution for ([T]). To formulate the 
result we use the following notations for slow variables: 

X j — X ^ t j — S"^ t ^ J — 1^2. 



2 



Let us define 
and 



K^2,h) = L[uot2 + kx2 + 5"]. 
The main result of the paper is formulated in the Theorem. 

Theorem 1 In the domain I < the formal asymptotic solution of (QP has 
the form 

U {x, t, e) ~ eui{xi, ti, ^2) exp{i{kx + ut)} + c.c. 
r /2 ] 

The amplitude Ui = ^'exp < —i — G(x2,t2) } , where G is an antiderivative 

[ Au J 

of 

1 1 



9[X2,t2) 



L[ujt2 + kx2 — S] L[ujt2 + kX2 + 5'] ' 

with respect to t2- 

Function \l/ is determined by the nonlinear Schrddinger equation 

ia;ai2*-9j^ + 37|*P* = 0, ( = uxi + kti. 

In the domain I > the formal asymptotic solution has the same form 

U (x, t, e) ~ evi{xi,ti,t2) exp{i{kx + ut)} + c.c. 
r /2 ] 

The amplitude Vi = \E'exp <i — G(x2,t2) } and is determined by the non- 

l 4a; J 

linear Schrddinger equation also and initial datum on the curve I = 
^{t2, C)b=o = e—i^ + — .'^ ^, 

2r(i-zV) 

where ^ = ^'(^2, C)l«=-o 

3 Numerical Simulation 

Here we illustrate the analytical results which are formulated in Theorem 
[H Let us consider ([T]) with 7 = —1/3 and S{t2) = t\/2. We consider the 
solution when /c = 0. It yields u = 1 and = xi. The envelope function is a 
solution of the nonlinear Schrodinger equation: 

2idt^m -dl^> - l^p^ = 0. 
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It is well-known that this equation has a soliton solution: 



COSh(?7(Xi + Kt2)) 

In this section we reduce the problem and take into account the envelope 
function as 77 = 1 and /t = 0. In this case the asymptotic solution before 
resonance has the following form: 

cosh(a;i) 

In this simplest case the resonant curves are ^2 = ~2, t2 = and ^2 = 2. 
Let us consider e = 0.1 then the resonant curves are t = —200, t = and 
t = 200. Following figures show the behavior of the solution for ([T]). Initial 
data are represented as: 



''^|i=-900 



Vcosh(£x) 

dtu\t=-goo = edt — —. — re ^ « vt^t^;; 
\cosh(exj 



-900 



-900 
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Fig.l. 3-Dimensional behaviour of numerical solution of ([T]) 



5 




Fig. 2. Profile of tlie numerical solution when x = 

4 The first external expansion 

In this section the formal asymptotic solution is constructed in the domain 
before the resonance. This domain is defined by the condition / < 0. The 
asymptotic expansion has the form of WKB-type. The leading-order term of 
asymptotic expansion has the order of e and oscillates. 

Theorem 2 When I < the formal asymptotic expansion for the solution 
of (Cy modulo 0{e^) has the form 

N 

U=Y.e^Un, iVGN (2) 

n>l 
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where 

Ui = ^ ni(j,o)(i^i,a;i,t2)exp{«j(A;a;2 + a;t2)A^} (3) 
i=±i 

and 

Un= Un,{j,m)iti,Xi,t2)exp{iXj,m{x2,t2)/e^). 
{j,m)£Un 

Here Xj,m = j{kx2 + ujt2) + mS is the phase of the ocsillating mode and the 
set fin contains the pairs (j, m) such that: 

fii= {(±1,0)}, 

fin = {(j„-i, m„_i ± 1), (j„_i, m„_i) G fi„_i} U 

{{Jh + Jh + Jh^^h + +^h)^ 
h + l2 + h = G ili^Jq G N}. 

r/ie coefficient Ui{i,o) defined out of [W^) . The higher-order terms Un^{j^m) 
are defined by linearized equation [T^) as {j,m) = (±1,0) and by algebraic 
equations p7| ) as {j,m) ^ (±1,0). 

This expansion is valid in the domains 

-E''L[Xj,rn] > 1, V(j, m) G uJLifi,. 

This theorem give the asymptotic solution for a set of the domains on the 
plane {x2,t2)- But this theorem does not give an answer about a connection 
between the asymptotic solutions for nearest-neighbour domains. 

This theorem does not give the solution in the form asymptotic series for 
([T]) in general. Because the resonant lines L[xj,m\ = are dense as (j, m) G 
and full series ([2]) does not asymptotical anywere. 

Let us consider an example as S* = ^2/2 and k = 0. The set fi„ contains 
a phases Xj,m as j = 2g ± 1 and q,m = ±0, 1, 2, . . ., \j\ + \m\ < n. In this 
case uj = 1 and the lines 

l-{2q + if ± m[-2(2g ± l)t2 - mtl] = 

are resonant. These lines 

are dance everywhere when N = 00. But for any segments of ([2]) the sets of 
resonant lines are finite. When = 3 the set is t2 = —2,^2 = and t2 = 2. 
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4.1 Construction of pre-resonant solution 
4.1.1 Derivation of equations 

Let us substitute ([2]) in equation ([T]) and gather the terms of same order of 
e. As a result we obtain a recurrent sequence of the equations: 

dlUr, - dlU^ + Un= -f COS )f/„_l - 

j+l+m=n 

-2dt,dtUn-2 - 2dtA2Un-3 - dlU^-^- (4) 

Let us consider this equation at n = 1: 

Lu'^ + k'^ + 0) exp{i(kx + cut)) + c.c. = 0. 

This equation defines the dispersion formula: 

,2 _ l2 , 1 



UJ 



Equation at n = 2 is 



[~ (^t2Xi,m)^ + {dx2Xj,m)'^ + lj«2,(j,m) exp (zXj.m/^^) 

(i)™')Sf22 

-^2ia;9tiUi,(i,o) - 2z/i;(9a;iUi,(i,o)^ exp(z(A;x + ut)) - 

f cos(S'/e^ + 0o)wi,(i,o) exp(i(/ca; + cut)) + c.c, (5) 

where ^2 = {j,fn), j = ±1, ''^ = ±1- 

We construct a bounded solution respect to fast time t. We remove terms 
with modes exp{i{kx + ut)) from the right-hand side of ([5]). It yields: 

t^5ti^ii,(i,o) - kdxiUi,(i,o) = 0. (6) 
Solution of this equation is an arbitrary function with respect to ( = kti + 

UXi. 

The set of of Q2 is defined by modes of right-hand side of The 
amplitudes for these modes are given by following formulas: 

( - {dt^Xj,mf + {dx2Xj,mf + l)^2,(j,m) = -^/«1,(1,0), (7) 

where (j, m) e {(1,±1)} and complex conjugated equations for (j, m) G 
{(-1,±1)}. Thus we define U2. 
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It is easy to see that U2 has singularity of first order on the curves: 

L[{kX2 + UJt2) + S] = 
L[{kX2 + Ujt2) - S]=0. 

Equation as looks like 

d^Us - dlUs + 1/3 = +2iujdt2Ui^(^i^o) exp{i{kx + ut)) 

+ 22^42 Xi,m<9tiM2,(i,m) " 

{j,m)en2 

'^id^2Xj,mdxiU2,{j,m) exp(?Xj,m/£^^) + 

(5*^1^1,(1,0) - 5^iMi,(i,o)) exp{i{kx + ut)) + 
fcos{S/e'^) J2 ^2,(j,m)exp(zXj,m/£^^) 



(9) 



+7 



ul ,^ Q-> exp{3i{kx + ujt)) + 



3|Ml,(l,0)rMl,{l,0)eXp(2(/cX + LUt)) 



+ c.c. = 0. 



(10) 



This equation has a bounded solution with respect to t if the right-hand 
side does not contain the terms with exp(±2(/cx + ut)). Thus we obtain two 
equations. First of them defines the dependence of M2,(±i,o) with respect to 
xi and ti. 

^dt,U2Xi,o) - kdxiU2,{i,o) = 0. (11) 
The second equation defines the dependence of Ui (1 0) with respect to ^2 and 
C: 



371^1,(1,0) r^i, (1,0) = -^9{x2,t2)ui,(ifi), 



where 



9{x2,t2) 
Let us denote 



+ 



L[ujt2 + kx2 — S] L[u!t2 + kX2 + 5*] 



wi(i,i) = exp I ~"^^<^(^2, ^2) I , 



'121 



where G is an antiderivative of g{x2, ^2) with respect to ^2- Then the function 
is determined by the nonlinear Schrodinger equation 



2ztu9t2^ - dc"^ + 371^1^^ = 0. 



3 = ±1- (13) 



The coefficients with exp(z(A;22; + ujt2 ± S)/e'^) and their complex conju- 
gations are defined by: 

The coefficients with exp(z(A;a;2 + ujt2 ± 2S)/e^) and their complex conju- 
gations are defined by: 

W3,(i,±2) = ~^^[~~^-^^2,(i,±l), i = ±1 (14) 

The following equations and their complex conjugations define the term 
with coefficient exp{3{kx2 + u!t2) / e"^) and exp{—3{kx2 + ujt2) / e"^) accordingly. 

1 

""3,(3,0) 



(15) 



As a result one obtains f/3 where 

n-s = {(1, ±m); (-1, ±m), m = 1, 2; (±3, 0)}. 

The third correction has the singularity of the second order with re- 
spect to L[{kx2 + uit) ± S] on the curves ([H]) and ([HI), and the singularity of 
the first order with respect to L[{kx2 + ^t) ± 25*] on the curves: 



L[{kx2 + ujt) + 2S] = 0, L[{kx2 + tot) - 2S] = 0. 



(16) 



Let us consider equation (j4j) for n > 3. The right-hand side of (j4]) contains 
the phases Xj,m, where (j, m) G r2„ and 

1]„ = Qn~i U {(j„,_i,m„_i ± 1), (j }u 
{{jii + ih +jh,'m-h +^i2 +mis)Ji + l2 + h = n, {ji^,mi^) G ili^Jg G N}. 
The amplitudes are defined by the following formula as (j, m) ^ (±1,0): 

_ 1 / , 



2^C^t2Xj,m^ti^n— l,(j,m) 2idx2\dx^Un—l,{i ,m) 



If If 

2 ~ ^"-1^0^™+!) 
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E 



^qi , (ji ,mi ) ^52 , (i2 ,m2) ^93 , (js ,"^3) 



(17) 



gi + 92 + gs = 

j = Jl + J2 + J3 

m = mi + m2 + ms 
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If (j, m) = (ihl,0) then the functions Un-i,(±i,o) = ^n-i,(±i,o)(^25 C)- Function 
^n,(±i,o)(^2 5 C) stay indefinite and Un-2,(±i,o) is a solution of the equation: 

2iujdt2Un-2,{±lfi) + <9^Mn-2,(±l,0) + 

7 ^ '^gi,(ii,mi)'W(32,(i2,m2)'^(j3.(i3.m3) = 0. (18) 

qi + q2 + q3 = n 
Up, rup) e 

jl +32+ js = ±1 
mi + 777,2 + ^3 = 

The term f/„ has the singularity on the curve „i] of the order p = 
ra — (|j| + |m|) + l. It means that the series for the solution loses the asymptotic 
property in a neighborhood of this curve. 

The series for the solution has an asymptotic property if 



« 1. 



It yields the domain of validity for the constructed asymptotic solution: 

^L[Xj,m] < 1. 

The Theorem |2] is proved. 



5 Internal expansion 

5.1 Expansion near primary resonant curve 

The primary singularity is situated at the curves = 0. It is easy to 

see that the ra-th coefficients of expansion ([21) has the singularity of the order 
n — 1 on these curves. 

A typical local resonance generates new harmonics with 

ki = k + 5^.2'S'|l[xi,±i]=o, as \ki\^k 

and 

Ui = U + <9t25'|L[xi,±i=0- 

That case may be studied by the same approach as [9]. 

Here we consider a special case of the local resonance when 
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This special case we call by the local parametric resonance. In this case a 
new harmonics are not generated in the leading-order term of the asymptotic 
expansion but the envelope function changes. 

In this subsection we construct the formal asymptotic solution near = 
and obtain an connection formulas for the nearest-neighbor solution ([2]). 

Theorem 3 In a neighbourhood of the parametric resonant curve the formal 
asymptotic expansion for the solution of (QP modulo 0{s^~^^) has the form 

N 

U{x,t,e) = Y,e''Wn{xi,ti,X2,t2,e), (19) 

n=l 

where 

" (2j — 1)5" 
Wn= 'Wn,j{xuti)exp{i — ), 

j=-n+l 

Function Wnj is determined from equations l[21\) . [23^) and l[2^) . This expan- 
sion is valid when 



5.2 Formal Construction 

In this section we prove the Theorem [31 
Substitute ([19]) into (P. Define: 



e\ = -\{d,,Sf + \{d^^Sr + l. 
Gather the therms of same order of e. In order of it yields: 



(20) 



W2,j expi^z- 



E 

i=-ijyo,i 
.(2j-l)5, 



2^2 



1 35* 1 3S 

-wi^ifexp (z— ) + -wi^^ifexp ( - i—) 

-wi,i]exp(z— ) + 

^uJT3T]exp(-z^), 
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This formula contains linear independent exponents with respect to fast vari- 
ables. Collect the coefficients of such exponents. As a result one obtains 
equations for wi,±i 



/ 



and a formula for W2 ±3: 



In order of we obtain the following equation: 



:2i] 



E 

j = —n + 1 
0, 1 



l((2j- 1)9^,5)2 + 1 



(2j - 1)^. 



1 s 

-/M);^]exp(i— ) + 

[{id.j,^Sd.j,^ - dt2Sda;^)Wn~l-l + \Wn-l-l + 

1 „ , , . s . 



i,_ijexp^-2— , 



t2 "^n-4,2j-l 



+ 



7 E 

p + g + r = n 
jp + Jg + Jr = 2j - 1 



E 

j = — n + 1 
0, 1 

2(2j - l)dt^Sdt2Wn-2,2j-l + 

/ / 

^'U^n.-l,2j-3 + -^Wn~l,2j+1 + 

.(2 j-l)g . 
2e2 ^■ 



(22) 



Collect coefficients with linear independent exponents with respect to fast 
variables. It yields: 



idt^Sdt^Wn,! - id^^Sd^^Wn,! + XWn,l + ^fWn,l = F^,! (23) 
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where 



Fn,l = -idt^Sdt^Wn-1,1 + idx^Sdx^Wn-l,! + 



(24) 

ni + n2 + = n + 1, 
ki + k2 + h = 1 
kj G J j — 1)2,3 

The term Wnj, j 7^ 1 is determined by algebraic equation 

{dt2SYwn-2,j - {dx2Sf^n~2,j " df^Wn~2,j + dl^Wn-2,j " 
2dt2dt^Wn-3,j + 2d^^d,^,Wn-3,j dlWn-4,j + dl^Wn-4,j - 

Yl 'Wn,M^n2M^n3,k3j ■ (25) 

'^l + ^^2 + ^^3 = + 1, 

ki + k2 + k3= j 
kj G ^Uji j — I5 2, 3 



5.2.1 Characteristic variables 

The function w„ i satisfies equation fl2T]) . The solution is constructed by the 
method of characteristics. Define the characteristic variables a, ^. We choose 
a point (x5,t?) such that C?x2^l(x0,t0) 7^ as an origin and denote by cr the 
variable along the characteristics for equation ( 12T]) . We suppose a = on 
the curve A = 0. The variable ^ mensurates the distance along the curve 
A = from the point (a;^,^?). This point (x?,t5) corresponds to = 0. The 
positive direction for parameter C, coincides with the positive direction of X2 
in the neighborhood of {x1,t^). 

The characteristic equations for (1211) have a form 

^ = dt^Siexi^eti), ^ = -d^^S{exi, eh). (26) 
da da 

The initial conditions for the equations are 

Xi|o-=o = 2;?, ti|o-=o = ^i- (27) 
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Lemma 1 The Cauchy problem for characteristics has a solutions when 
lo"! < Cie~^, Ci = const > 0. 

Proof. The Cauchy problem fl2U|) . fl?r|) is equivalent to the system of the 
integral equations 

ti=t1+ / dt^S{exi,eti)dC, xi = Xi— / dx2S{exi, eti)dC. (28) 
JO Jo 

Substituting ^2 = (^i — ^i)^^, 5^2 = (a^i — x1)e, we obtain 



to 



dt2 S {x2 - exi, ^2 - £t°i)dC, X2 



8x28 {x2 — £x1, t2 — £t1)dC. 



The integrands are smooth and bounded functions on the plane X2, ^2- There 
exists the constant Ci = const > such that the integral operator is a 
contraction operator when e\a\ < Ci. Lemma [T] is proved. 

It is convenient to use the following asymptotic formulas for the change 
of variables (o", 0- 

Lemma 2 In the domain |cr| ^ the asymptotics as e ^ of the solutions 
for Cauchy problem / f^) . [2l\ ) have the form 



N 



- xliO = -adx,S + ^ e"a-+'gn{exr,etr) + 0(£^+V^+2)^ 

n=l 

N 



(29) 
(30) 



n=l 



where 



da"" 



(7=0 



da'^ 



a=0 



The lemma proves by integration by parts of equations ( l28i) . 
The next proposition gives us the asymptotic formula which relates the 
variables a and A as a, A — ±00. 

Lemma 3 Let be a <^ e~^, then: 

A = (^(0^ + 0(6^2), a ^00, ^{0 = ^ ■ 

da cr=o 

Proof. From formula (l20l) we obtain the representation in the form 

00 
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where 



It yields 



Xj{xi,ti,e) 



1 



A(xi, ti, e)\a-=o- 



Let be 



d^l 



.d^X. 



da^ 



> const, ^ G -R. 



The function dX/da is not equal to zero 
dX 1 



da 2 

Let us suppose dX/dcr > 0. It yields 



-d,M.2S + dt,Xdt,S] ^0 



(7=0 



The lemma is proved. 

5.2.2 Ordinary differential equations for local parametric reso- 
nance 

The equation fl^Tl) may be considered as an ordinary differential equation 
along the characteristic. 



I— h X{a,^,e)wi^i + -wi^i = 0. 

da 2 



(31) 



Let us expand A as a series with respect to e and consider (13111 as a pertur- 
bation of: 

dwi 1 /_ 

+ </?(4)o"^f^i,i + 0^1,1 = 0- 
aa I 

Let us change variable: k = aa, F = //(2a) where Q!^(C) = ^{0, then we 

.dwi 1 



obtain: 



dn 



+ KWi^i + Fwi^i = 0. 



(32) 



The same equation for Wn.i has the following form: 



. dWn,l 

dK 



1 



HWn,l + FWn,l = -Fn,l + - V Ai+i(-) Wj^i- (33) 

j+t=n— 1 
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A solution for the equation flH^ has obtained in [TU]: 

V{k- C) = ^(-1 - if-'^ FCD,{e-fV2K) + 



2 v27r 



1 _ + 2(-l - z)-i-VeV ^ C (e^V2, 



Fr(-2f 



where D^iy) is a parabohc cyhnder function, 2; = -i^^ — 1 and C is an arbitrary 
complex constant. 
Let us define 

wi,i = V{K,C). (34) 
Let us consider the nonhomogeneous equation: 

i'^ + kU + FU = G{k). (35) 
an 

To solve this equation we use two linear independent solutions of homoge- 
neous equation. These solutions are: 

Vi = d^V{K,C), V2 = dpV{K,C), C = a + if3. 

Let us define the Wronskian of these solutions as: 

W = Vi% - V2Y1. 

It easy to see that W = const 7^ 0. 
The solution of (1351) has a form: 



U = V, J\gV2-V2G){x)^ + 

V2 f{GVi-ViG){x)^ + ciVi + C2V2, (36) 

where ci and C2 are real constants which are parameters of the solution. 

The general formula (l36l) allows to us to solve the equations for any high- 
order term. 



5.3 Asymptotics as A ^ 00 and domain of validity of 
the internal expansion 

The domain of validity of the internal expansion is determined by the asymp- 
totic behaviour of higher-order terms. In this section we show, that the n— th 
order term of the asymptotic solution grows as A"~^ when A — 00. This 
growth of higher-order terms allows us to determine the domain of validity 
for internal asymptotic expansion ( fT9l) as X 00. 
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5.3.1 Asymptotic behaviour of leading-order term 

The leading-order term has the following behaviour as /t — > — oo: 

^i,i = Ci,ie<^-^i'^(-'^))(l+0(«:-i)), Ci,ieC (37) 

and 



2 2 
a TT ■ a 



V 11 — t^) / 



r(i-^f) 



as K ^ oo. 

This formulas were obtained using the asymptotic behavior of solution 



5.3.2 Asymptotic behavior of higher-order terms 

The n-th order term has the following asymptotic formula: 

Then the internal expansion is valid as |k| <^ or the same as £:|A| ^ 1. 

6 Matching 

The domains of validity for the internal and external asymptotic expansions 
are intersected. It allows us to match these expansions and as a result to 
obtain an uniform asymptotic expansion that is valid in both domains. 

To match the expansions one should reexpand the external asymptotic 
expansion in the terms of internal variables ^ and k, and equate the terms 
with the same order of e. It yields the indefinite functions in the internal 
asymptotic expansion. Here we make this step for the main terms of the 
internal and external asymptotic expansions before the local resonant layer: 



e[ k) exp{iS/e ) + c.c. 

P 



e[^!{t2X)Q^v{-i^G{x2M))Q^v{i{kx2 + ^t2)/e^) + c.c^, (39) 

as 1 < -K < £-1, {t2,X2) e L[xi,i] = and C = tux?(0 + kt\{i). This 
formula define the function wi-i{S,, k) as n ^ —oo. 
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The same formula for the matching is vahd after the resonant layer: 

\Ef(t2,C)exp(-2— G'(x2,t2))exp(z(/cx2 + ujt2)/e^) + c.c. j ~ 

e (^wi _i k) exp(2S'/5^) + c.c.^ . (40) 

as 1 < K < (^2,2:2) G = and C = u}x\{i) + kt\{i) . This formula 

defines the function \E'(t2,C) coming out of the resonant layer. 
The matching gives a jump for the ^ at the resonant line: 



^^(^2, C)U[xi,i]+o = e— ^'(t2, C) |l[xi.i]-o + 



1 +i)e^e^fe'"(')/v/? 



-^(t2,C)|L[x.a]-o) (41) 



6.1 Post-resonant expansion 

The asymptotic structure of the post-resonant solution has the similar struc- 
tire as the pre-resonant asymptotic solution. The main term of the solution 
is: 

U (x, t, e) ~ i(xi, ti, ^2) exp{2(/cx + ut)} 
{ P 1 

The amplitude Vi = exp < i— G(x2, ^2) ^ and ^ is determined by nonlinear 
Schrodinger equation also and initial datum on the curve = 

, (l + 2)e^e^«''^(2)jy^. 



^(^2,C)U[xi,i]+0 = e— ^(t2,C)U[x,,i]-0 + ^ _ •£/ ^fe,C)|L[xi.x]-0- 
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